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Abstract 

We introduce a technique for recovering a sufficiently smooth function from its ray transform 
over a wide class of curves in a general region of Euclidean space. The method is based on a 
complexification of the underlying vector fields defining the initial transport and recasting the 
problem in terms of complex- analytic function theory. Explicit inversion formulae are then given 
in a unified form. The method is then used to give inversion formulae for the attenuated ray 
transform. 

Keywords: Explicit inversion, complex analysis, transport equation, quasiconformal, harmonic cal- 
culus, attenuated ray transform, Riemann-Hilbert problem 

1 Introduction 

In several engineering situations one deals with data consisting of the line integral of a function 
and the goal is often to recover that source function from its integral over a class of lines. In the 
arena of medical imaging, this arises in positron emission tomography (PET), single photon emission 
tomography (SPECT), and (originally) CT-scan tomography [9j. In other applications (in geophysics 
[29J and non-destructive electrical imaging techniques [5] such as electrical impedence tomography, 
EIT) the line integral is instead taken over a class of one-dimensional curves in either Euclidean space 
or more generally, a Riemannian manifold. This type of data is referred to as a ray transform. 
In geophysics for instance, the problem can arise as the linearization of determining geophysical 
properties of the Earth from travel-time measurements [29]. Quite often the physics will also dictate 
that the signal suffers some absorption along its trajectory and is attenuated, the data then called, 
not surprisingly, the attenuated ray transform. 

The mathematical applications, properties, and uses of these integral transforms and their inverses 
are discussed in great detail in [TJ [131 O 122] and include harmonic analysis, algebraic curves, tensor 
geometry, and partial differential equations to name a few. Generally, explicit inversion formulae 
over curves other than lines (geodesies of a Riemannian manifold |18| . say) tend to restrict focus 
to manifolds with a strong amount of symmetry [HI [TH [I5j [U [27] and do not include the effects of 
absorption encountered during propagation. An exception to this statement can be found in |19t I21j. 
We restrict our attention in this paper to curves in a 2-dimensional region of space. 

The method used in this paper generalizes a technique first used in |23] for lines in Euclidean 
space and later generalized in [3] for geodesic rays in hyperbolic geometry for giving an explicit 
inversion formula for the attenuated ray transform in each case. The technique we present rests on 
the complexification of a certain class of differential operators in IR 2 which allows us to recast the 
problem in terms of complex analysis in the unit disc. Once the problem is cast in this light, we 
use the classical Poisson formula |20| to give us a reconstruction formula. Excellent introductions to 
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complex analysis and conformal mappings are [251 El ED] and the classic [1] . Good introductions to 
quasiconformal mappings and Beltrami equations (and their generalizations) can be found in [2"1 124] . 
References on Blaschke products and multivalent mappings can be found in [61 [10] . 

An outline of the paper is as follows. In section [2j the general setup, notation, and a quick review 
of the essential operators used throughout the paper are presented, together with the main result we 
are seeking to establish. In section [3] we begin the complexification procedure by introducing a new 
(complex) parameter A into the transport equation introduced in section [2] and give a classification of 
the vector fields under consideration as those of type H. Much of the heavy lifting is done in the more 
technical section [4] where we find and analyze the Green's function of the new parameterized complex 
partial differential transport equation. We will establish that condition H is sufficient to guarantee 
holomorphicity of the solution of this equation in terms of the new parameter A. We evaluate the 
asymptotics of the solution as our complex parameter A tends to the unit circle from both inside and 
outside, i.e. as |A| — > 1 T and see that in fact its imaginary part depends on the data we are interested 
in. Once this is established, we use this fact in section [5] to give our desired reconstruction formula in 
the non-attenuated case. The rest of section [5] uses the non-attenuated formula to give an integrating 
factor solution for the attenuated case, which requires an additional constraint to condition H. We 
offer brief concluding remarks in section [6l 



2 Preliminaries 

The Stationary Transport Equation 

We let 7 : M? 3 (t,s) *— > ~f(t,s) £ C I 2 be a real-analytic diffeomorphism where CI is an open, 
bounded, simply-connected region of the plane (a domain). We consider 1R 2 = C by the standard 
isomorphism so that 7 is identified with 7 1 (i, s)+ij 2 (t, s). Then, (w, w) are (independent) coordinates 
on CI where w = ~f(t, s). Because 7 is a diffeomorphism, its differential is injective and therefore induces 
a vector field on Cl via its differential under the rule (cj)*X)(f) = X((p*f). Consider 7*Jr. This gives 
a non-degenerate field of the following type; 



d _ d 

X\ w = [jl(w)- \- n(w)— w£ CI, \fi\ > 

ow ow 

which acts on pushforwards in w of functions on Cl and the non-degeneracy is ensured by the regularity 
of the curves -y(t,s). The equation of interest is the stationary transport boundary value problem 
X\ w u(w) = f(w), for w 6 Cl, f(w) € Cq°(CI) with lim^-oo u(w(t, s)) = 0, i.e. the BVP 

u\ d _n = (2) 

Anticipating our desire to complexify the above equation, we be wanting to exploit the 50(2) sym- 
metry of the unit disc which is a priori not available to us in this more general domain by appealing 
to the Riemann mapping theorem [11]. Denote the unit disc by D + = {zsC:[z|<l}, the unit 
circle by T = {z € C : \z\ = 1}, and D~ = C/{D + U T}. Let K(z, w) be the Bergman function of the 
domain CI, then the function 



Z{W) = j c \jKjtl) K{tX)dt (3) 

gives the unique biholomorphism mapping 0, into D + , the unit disc, with z(() = 0, -z'(C) > as in 
|22j and (t,s) give coordinates on D + through composition since 7*2; maps M 2 into D + . Because of 
this equivalence between our initial domain Q and the unit disc all further results will be presented in 
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the disc. If £1 was all of M? (and the Riemann map was consequently unavailable) the method below 
will still work since R 2 has the needed rotational symmetry. 

We therefore use (z, z) as coordinates on D + and have a new vector field on D + given by X\ z = 
z*X\ z r w y and fx — > o z~ 1 and likewise for p,. By slight abuse of notation we denote {z*[i}-§^ o 

z^ 1 by n(z) and {z*fi}§= o z~ l by fj,(z) so that field of interest is 

X\ z = Li(z)—+fi(z) — , z£D + 

We define t{z) = z^w^t and s(z) = z^w^s, smooth functions on D + . 

The method of characteristics gives the following solution to the transport equation X\ z u{z) = 
f(z), u(z(— oo,s)) = as 

u(z) = {Dxf)(z) = \ ! f(z(t , s))sign(t(z) - t Q )dt Q (4) 



and since the integral curves of X\ z are just the image of integral curves, i.e. y*z* = (z o -y)*. we 
define the ray transform of a source function f(z) over the integral curves of X\ z indexed by s to be 

(//)(*)= / f(z(t,s))dt (5) 

We will later be using the following extensions of these operators given below: 
Symmetrized beam transform 

(D g i>)(z) = \ f iP(e ie z(t , s(ze- ie )))sign(t(ze- ie ) - t )dt ^ G L l (D+) 

Ray transform 

m(s,e td ) = (l e ip)(s) = [ iP(e w z(t,s))dt V G L\D+) 



We will always use 9 and e interchangeably, its meaning clear from context. 

We will have occasion to use the Hilbert transform H of a function defined (see e.g. |30| ) as 
the following Calderon-Zygmund principal value integral operator 

(Hip)(x) = -p.v. [ ^-dy VeLP(R), P>1 (6) 
7T jRX-y 

1 — \z\ 2 

Lastly, we will be using the standard Poisson kernel of the unit disc given by P(z, 6) = 1 _ e -;9 2 |2 ■ 
We recall (see, e.g. [SET]) that the Poisson kernel generates harmonic solutions v(z) of the BVP 

Av = z G D + 

v\t = 9 

given by v(z) = ^ f T P(z,9)g(e t9 )d9. The Poisson kernel is also deeply connected to the study of 
inner functions c.f. [1U|, 125]. 

The main purpose of this article will be to show that given suitable conditions on /j,(z, z) and s(z, z) 
that 2 

f(z) = ^fj p (k, 9)X e A -H(I f)(s(ze- l9 ),e ie ))d9 i = 1, n 

where the \i{z) are functions to be introduced later. With this established, the above formula is 
used to give an integrating factor method to find a similar reconstruction formula for the attenuated 
ray transform along the same lines. The above is a type of inversion formula known as a filtered 
backprojection type [9]. The procedure used to derive the above main result can be best thought of 
in the following heuristic scheme 
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1. Model: Writing down the linear stationary transport equation for the dynamics 

2. Symmetrizing: Introducing a rotation parameter A = e td into the integral curves of the trans- 
port PDE 

3. Symmetry-Breaking: Complexifiying the parameter introduced in step 2 by moving A "off- 
shell", i.e. |A| ^ 1 

4. Analysis and Asymptotics: Evaluating the dependence of solutions to the complexified equa- 
tion on our parameter A and examining limiting behavior 

5. Reconstruction: Using holomorphicity of the solutions to write the inversion formulae as 
Poisson integrals of their asymptotic boundary values found in step 4 

The reader may find some benefit from keeping the above rough outline in mind throughout the 
following. In this section, we have finished step 1. Steps 2 and 3 are handled in the next section. Step 
4 is done in the more technical section 01 and the final step is given in section 

3 Complexification of the Transport Equation 

Since D + is acted on transitively by 50(2) we will define the conformal map A : (z, z) — > (Xz, jz), for 
A G T the unit circle. Notice that if $(-,s) is a set of integral curves of D + , that z (A*$(-, s)) are 
conformally related curves in f2. Then, for A G {D + U D~}/{0, oo} we consider A*-X"| 2 = X\ to be the 
so-called "complexification" of X\ z . We remark that \*X\ Z takes the form /x(f , Az)A^ + ^)\~§i 
or 

X x = Z(z,\)— + p(z,\)— XeD ± /{0,oo} (7) 
oz oz 

with j£(z, A) = p(z, A) = X*fj,(z) and Xp(z, A) = pl(z, A) = A*/x(z). We also define X^ = ±i(— A) J| + 
p(z, A)JL) as a vector field orthogonal to X\ when A = e ld . Namely, Xg ■ Xq = ±(£(z, e ld ), p(z, e td )) ■ 
(-i£,(z,e ie ),ip(z,e ie )) = ±i(\£(z, e id )\ 2 - \p(z, e iB )\ 2 ) = in the standard inner product (•,•) : C 2 -> C. 
The factor of i is needed to make X^u real- valued and the choice of ± is determined by whichever 
satisfies the condition X^s > 0. Since X^~ = a(z)z*J^ for some real-valued a(z), this determines X± 
uniquely. Since we could just as well reparameterize with —s we will, without any loss of generality, 
avoid keeping track of signs by just assuming that X^ = i(—£{z, A)J| + p(z, A) J=). 

We likewise define s(z, A) and t(z, A) as X*s(z) and X*t(z) respectively for A G D =t /{0, oo}. A word 
on notation: ^| and k z are equivalent, as are ^= and k z , and we will use them interchangably. 

We remark that equation (|7|) has no direct physical meaning since the complex parameter A, when 
taken to lie away from T = dD + , is in some sense artificial and may be best thought of as a complex 
parameter indexing a class of complex partial differential equations given in ([7]). 

Next we reduce the scope of our consideration to the class of vector fields X\ so constructed to 
consist only of those of type H: 

Definition A complexified vector field X\ = a(z, A)Jg + b(z, A)J=, induced in the manner above as 
A*X| Z , A € -D^/jO, oo} from a real field X\ z , is said to be of type H if the following holds: 

• a(z, •) is a holomorphic function of A for A G D + and has at least one zero A = Xi(z) G D + 

• b(z, •) is a meromorphic function of A for A G D + and has no zeroes in D+ 

• b(z'-) ^ S a holomorphic function of A for AgD + and has at least one zero A = Xi(z) G D + 

• s(z, •), ; d s 0v) are m eromorphic functions of A for A G 
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where, as in the above, s(z,X) = A*s(z) is the complexified transverse foliation parameter of the 
integral curves of X\. 

We are, in the above, treating z and A as independent variables, and holomorphicity is to be 
thought of in the standard way of functions of several complex variables |16j . We are therefore not 
requiring any of the above functions to be holomorphic in the z variable. 

Because ^Jl'jI * s holomorphic in A € D + its zeroes are isolated. Also, since * s holomorphic for 
A € D + and since conformal mappings map boundaries of Jordan domains into boundaries of Jordan 



domains (see |22]), then _^' e , J = ^£4 for some y G T and thus I ^ z '^\ 



= 1. Since we assumed 

|A|=1 



that there is at least one zero Aj, the maximum principle implies that | ffi^ [ < 1 for A € D + . We 
then get the following simple 

Lemma 3.1 ^Jf z ' } \) ^ as a fi n ^ e number of zeros, \i(z) with multiplicities rrii(z) 
Proof This is a simple consequence of the argument principle (fllf)- Namely, one has 

d t(z,\) 



and \\ ^ s holomorphic, hence so is ^ p \z x) > on ^ e re 9^ on D + . They are also both continuous on 

dX p(z,X) I 



T. Therefore, |£%r}| < M < oo for A € D+. Thus, 



2- 



? m - £ i | l, rf l!f' ,(2 ' A)dA|< iX MdS=M (8) 

Henceforth Xi will always be used to indicate a value in the unit disc for which ^f^y (and £) vanishes. 
The bounded holomorphic functions mapping the unit disc onto itself and having a finite number of 
zeroes can be uniquely written as a finite Blashke product (c.f. [6j [10] ) so that p^'x) can be given in 
the form C(z)n? =1 ( A ~ A ' ) m ' with \C(z)\ = 1, and with nii and Aj possibly depending on z. 

Furthermore, since l^f^yl < 1 f° r A € D + we also have that the complexified transport equation 
X\u(z, A) = f{z) can be rewritten as 

u- z (z,X) = ^f^u z (z,X) + ^- (9) 

which is a forced Beltrami equation as in |17t H] . 

4 Solving the Complexified Equation 

In trying to solve the complexified transport equation 

X x u(z,X) = f(z) (10) 

we will again be changing variables. Notice that X\s(z, A) = on that region so that s is still a constant 
of the dynamics. This is obvious from the fact that integral curves are mapped by diffeomorphisms 
to integral curves, however to be precise, when |A| ^ 0, 

d d 
X\s(z,X) = A*X| z A*s(,z) = A*z*w*— X^z^WtS = (A o z o w )*q^ (^ ° z ° w )* s 

ds 

= (X o z o w)*— = 
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since s and t are independent coordinates. Thus 



..ds(z, A) , .^ds(z,X) 
g(z,A) V ' +p(z,X)- 







(11) 



<9z dz 

The Riemann removable singularities theorem [11] applies when A = 0. We will need the following 

Lemma 4.1 On < |A| < 1 the Jacobian ds(z) = \s z (z, A)| 2 — |sg(z, A)| 2 is positive 

Proof Since (t,s) i— > (z, z) is a diffeomorphism and A : (z, z) i— )■ (-|,zA) is conformal on < |A| < 1 



&s(z,A) dt(z,X) 
8{zX) d{zX) 



(12) 



so f/iaf 



.ds(z,X)dt(z,X) ds(z,X)dt(z,X). 



dz dz dz 
implies |s z (z,A)| 2 ^ 0. Then, 

ds(z) = \s z (z,\)\' 



< \s z (z,X)\(\t- z (z,X)\ + 



p(z,X) 



Xjz,X) 
p(z,X) 



s z (z,X)\ 2 > \s z (z,X)\ 2 (l 



, gz,X) 
p(z,X) 



Uz,X)\) 



> 



since < 1 for X € D+ . 



Since X\s(z, A) = 0, s*ATa = s*X\s(z, A) J|. We are interested in solving X\G\(z; Zq) = 5{z — Zq) 
and we can achieve this by solving s*X\s(z, X)-g=(s*G\) = \ds(z)\5(s(z, X) — sq). We therefore need 
to compute the term s*X\s(z, A). To this end, 



implies 



whence 



Oz 



az,x) = - P (z,x) 



dz 



ds(z,X) 

ds(z,X) 
dz 



ds(z, A) 
dz 



+ p(z,X) 



ds(z, A) 
dz 



ds(z,X) Q z 
dz 



dz 



ds(z,X) 
dz 
1 

Q(Z,X) 



(\s z (z,X)\ 2 -\s- z (z,X)f 



ds(z) 



with Q(z, A) = 



ds(z,\) 
-dz_ 



By recalling that |-| > 1 for |A| > 1 and going through the preceding lemma 



p(z,X) •> " i pi 

mutatis mutandis we see that ds(z) is likewise negative on D~ and hence the Jacobian of s(z, A) 
switches sign when A € so that our fundamental equation then becomes 



1 



d 



which gives G\(z; zq) 



Q(z, A) ds 
signjl— \X\)Q{zp,X) 



tt(s(z)-s(z { ))) 



s*Ga = sign{\ — \X\)S(s(z, A) — s(z , A)) 
or rather 



G\(z;z ) 



1 ds(z,X) | 



dz \ z 



tt(s(z,X) - s(z ,X)) 



X € D ± /{0,oo} 



(13) 
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so that u(z, A) = f D+ G\(z; Zq) f (zo)dfj,(zo) solves X\u(z, A) = f(z) for A G D ± /{0, oo}. We have used 
the fact that = S(z) as shown in 



Remark We will only make use of results in our formula which follow from condition H and thus 
results like (|13p are only used when A G D + . We will however present many results for A G D~ with 
the understanding that given an appropriate generalization of condition H (involving constraint on 
£ and p for A G C/D + ) the results are true. The advantage to this approach is it makes use of the 
symmetries and parallels of several of the formulae for A G D^. Thus, in the "— " versions of several 
results, condition H is necessary but not sufficient. 

With the above in mind, since 



d(t,s] 



dz 



d(z,z] 



ds 
dz 



dz ds dt 1 
didEdz^ 

az 



dz dt 
dt dz 



dz dt dz dt 

'^did^ + dt ~dT 
dt 

~ z *di 



we can rewrite G\(z; zq) as 



Then for ip G C^(D + ) and dfi(z) 



G\(z,z ) = -A 



d(t,s) 




d(z,z) 




n(s(z) — 


s(*b)) 



dzdz 
2i 



dxdy, the standard Lebesgue measure on 



(14) 

(15) 
C we have 



d(t,s) 

J\(D+) n{s(z) - 



Z() 



d(t,s) 




d(z,z) 





dp{zo) = / (A»— — — — z -f—-dfj,(z ) 
Jd+ tt{s{z) - s{z )) 

Iff \*ip(z(t ,so))dt ds 

2ni J R J R s- s 



(16) 



so that J D+ ip{zo)G\{z; Zo)dpi(zo) stays bounded since A fixes the unit disc. A similar argument works 
for A G D~ /{oo}. Because of the meromorphy assumptions stated in condition H, we have that when 
z / Zq, G\{z;zq) is a holomorphic function for A G D /{0, oo}. Since s(z, A) and s z (z,X) have the 
same order of (possible) pole at zero, G\(z; Zq) stays bounded even at A = and we get the following 

Proposition 4.2 u(z, A) is holomorphic for A G 

A similar argument applied to -§;G\(z; zq) and -§=G\(z;zq) shows that u z (z, A) and Uz(z, A) re- 
spectively are also complex-analytic for A G D^ 1 , a fact we will make use of in our final reconstruction 
formulae. 



Boundary Behavior 

We will be using the boundary values u(z, X)\\^t to arrive at a reconstruction formula. Therefore, 
ignoring the signum for the moment and letting ip G C Q yD (D + ) be a test function, then, using (|15p . the 
two-form ifi(zo)G\(z; Zo)d[j,(zo) equals 

-\*{ifj(\*z(t , s )) 1 rdt ds } 

2tti{s — So) 

so that we get the following 
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Proposition 4.3 u±(z, e s ) = lim D ± 3X ^ ei e u(z, A) = ^^(HIgf){s(e t9 z),6) + (D e f)(z) where the 
Hilbert transform H is taken with respect to the first variable. 

Proof First we examine s ^ z x)- s (zo A) w ^ en A = 1 — e (e << 1) and use the fact that s(z, 1 — e) = 
s(z, 1) — es(z, 1) + o(e 2 ) together with X\s(z, A) = to get 



Oil): 
0(e): 



and 



(^,1)|+^1)|)^1)=0 

(i(z, 1) A + p{z , l)^)s'(z, 1) = 1) A + ^ 1)^M*, 1) 



1)^ + l)§z)s(z, 1) = -(£'(*, 1) - p'(z, l)^-) Sz ( z , 1) 



-Z(z, i)s z (z,iy 



v 9A p^ 



A=l 



£(*>!) 



(17) 



so i/iat 



*im'(z,1) = -i£,{z,l)s z {z,l)- 



Ah 

OA p/ 



A=l 



5y a similar argument one can show 



Xiis'(z, 1) = l)s f (z, !)■ 



V <9A p/ 



A=l 



so i/iai 



Xiis'(z,l) = - 



(9_L) 




\ <9A pi 


A=l 



2 gfajj 



1) 



(18) 



Since | is given as a finite Blashke product C,{ z )^i=i{ \z\\^ ) mi ^ z \ we 



so that 



( 4-^1 



g t(*,A) 

see t/ifli 8 * ( ff 



E*>o m . 



l-IA.f 



j>0 J (A-A J )(l-A i A) 



and therefore 
Then we look at 



> 0, which, when combined with X^~s(z, 1) > giwes /rom jTgj) i/ia£ 

Xvs'iz, 1) > 

sign(is'(z, 1) - is'(z , 1)) = sign(t(z, 1) - f(2o, 1)). 

^(^(iojSo)) 



/ G 1 ^ e (z;z )ijj(zo)diJ,(z ) ^- -— / — — - 
Jd+ 27TZ J R J R s{z, 1 



e) - s(zq, l-e] 



dsodto 
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which as we have shown is 

1 f f ip(z(to,so)) 



2vri J R J R s(z, 1) - s(zq, 1) - e(s'(z, 1) - s'(z , 1)) 
Iff ip{z(t ,s )) 



dsodto 



dsodtQ 



+ o I I ~ s{z ,l))sign(is'(z,l) - is'(zo,l))if)(z(t ,so))dt dso 

1 JR JR 

1 H(I il>)(8(z),l) + 1 I sign{t{z,l)-t )i;(z(t ,s{z,l)))dt (19) 



2vri J R J R s(z, 1) - s{z ,l) 
1 
2 
-1 

2*"-*-'"-'^--"-' 2. 
XTuts, we /iaue 

u + (z, 1) = ^rH{I e ^){s(z), 1) + (L>i^) (20) 
For the general case, G e ze(z;z ) = G\(e~ %e z; e~ ld zq) shows 

u + (z,e* e ) = =±H(I e ij)(s(ze- ie ),e w ) + (D e ^(z). (21) 



An identical argument for U-(z, e ld ) shows that 



u±(z, e ie ) = T^H(I e i;)(s(ze~ id ), e w ) + (D e ^)(z) (22) 



5 Inversion Formulae 

5.1 No Attenuation 

We can now prove our main result. 

Theorem 5.1 If X\ is a vector field of type H, jfcQ 
then 

l-2n 



= fori = l,...,n and f(z) G C^(D + ), 



1 C 

f(z) = — J P(X u 6)X i tH(Ief)(s(ze- w ),e w )de 



gives an exact reconstruction formula for the density f based on the data I$f of ray transforms of f 
over the integral curves of Xq. 

Proof With P(z,9) = n~=^p > the Poisson kernel of the unit disc, and Cauchy's formula for holo- 
morphic functions, one has (120]) that 

i f 2n 

X x .u(z, ^) = ~r P(Xi,0)X e H(I e f)(s(ze- id ), e ie )d6 



4tt 



o 



1 f 2n 

+ — P(\i,e)X (D e f)(z)d8 (23) 
2tt Jo 



so that 

r-2iT 



X x .u{z,\i) = f(z) + — J P(X i ,e)X e H(I e f)(s(ze- ie ),e ie )de 



whereas 



X^u(z,Xi) = ^ J 2 J P{\ h 9)XtH{Ief){s{ze- ie ),e w )d6 

+ ^ fj P(Xi, 6)Xt(D e f)(z)d6 (24) 

Then since X\ = £{z, A) J| + p(z, A) X^ = A)J| + /9(,2, A) J|) and £(z, Aj) = ; we have that 

iX x .u(z, Aj) = X^ufo Aj) 
so that, on equating real and imaginary parts of t23\) and \2J$ , we get 

2vr 



n P(Xi,e)X^(D e f)(z)d9 = W P(X i ,e)X e H(Ief)(s(ze- ie ),e id )de 

and 

f(z) = J 27 " P(\i,9)X^H(I 9 f)(s(ze~ i6 ),e i6 )d8 (25) 



It's clear that formula (|25p could just as well be written in terms of the jump function (from the 
viewpoint of D^) 

<f>(z, e M ) = u+(z, e ie ) - u-(z, e M ) = iH(I e f)(s(z e - w ), e w ) 

as 

f{z) = —J P(X t ,9)X^-i<j ) (z,e ie ))de (26) 

an observation which will be useful in the next section. Recalling our previous remark about using 
only results from D + we could just as well use 

(j)(z,e ie ) = 2i%u + {z,e ie )) 

and remember that invoking D~ is only a useful mnemonic. 

5.2 Attenuated Ray Transform and Inversion Formulae 

We add a real-valued attenuation term a(z) € Co°(-D + ) to the complexified stationary transport 
equation to get 

(X x + a(z))u(z,X)=f(z) XeD ± (27) 
Using our Green's function G\{z;z$), we define 

h(z,X) = / G\(z; z )a(z )dfi(zo) (28) 
Jd+ 

and we use an integrating factor approach as follows 

e h(z ' X) X x u{z,X) + e hiz ' X) a{z)u{z,X) = e h(z ' X) f{z) 

so that 

X x e h ^u{z,X) =e h(z ' x) f(z) 
10 



whence 

u(z,\)= [ G x (z;z )e h ^- h(z ^f(z )d^z ) (29) 

JD+ 

Now, since 

h±(z,e ie ) = T^{HIea){s{ze- ie ),0)) + (D e a)(z) 

as before, we have the solution of the attenuated transport equation admits the following boundary 
values as |A| — >■ 1 T 

— h±(z e iB ) 

u±(z,e w ) = Te 2 , ' [g/4e fe± (-' e ' 9 )/}(.(z e ^),g)T^(ZJ e e fe± ^ e ' fl )/)(z)] 

— h+(z e ie ) 

= ^ ll [ gft{e #(W(^) /( yM0 }(g(ze -tf )|9) 
T ^(A^™^ 19 -^) /(-)e^ a )(-))(z)] 

Defining 

(J a , e /)( S ) = / e (/(-)e^°)»)( S ) (30) 
and recalling that Ig involves integration in t, not s (as does Dg) and therefore 

—h±(z e iB ) 

u±(z,e ie ) = T€ ' g(elr( g ( J ^)('(«- ifl 0,g)j a>o /)( a ( ze -<g) > g) + e -(^«)W(i3 9 /(.) e (^a)(-))( z ) 

_■ / 



so that 



(z,e ie ) = (u + -u_)(z,e w ) = - & ^ } g(e^ H ^ a )^ e " 6 '-)' e )/ a , e /)( g (ze-^), g) 



e 



2z 

e -(£> 9 o)W 



i7( e -ii^(^«)(-( e - ie -)^)/ aie/)(s(ze ^ ) ^ ) 

| e ^(/ s a)( S ( 2 e- 9 ),9) F ( e iff(; 9 a)( S ( e - w -),0) /([() ^ 



2i 

+ e -^H(^a)W.e-«),e) ff(e -iH(7 fl a)( a (e-«.),*)/ aie/ )}( a ( ze -ie) je ) 

We define C = cos( H(/ ^ ( f e ~^ e) ) and 5 = sin( g(J » a) (f e ~"^ ). Then 
e -(A>a)(*) 

e*) = [(C - iS)H{(C - iS)Ia, e f} + {C + iS)H{(C + iS)I a> ef}](s(ze- te ),6) 

= ie~^ a ^{{C ~ iS)H[(C - iS)I a>e fMze-* e ),8)} 

= ie^ D ^ z \CH{CI afi f){s{ze- i6 ),e) + SH{SI afi f){s{ze^ e ),e)) (31) 
= ie^ D ^ z \H a I afi f){s{ze-^),6) (32) 

where H a : f i-> CH(CI a fif)(s(ze^' l6 ),6) + SH(SIafif){s{ze" l6 ),6). We then can proceed in a manner 
similar to before since we have that e h<yZ ' X ^u(z, A) (along with its derivatives) is holomorphic and solves 
X\e h( - z ' X ^u(z, A) = e h ( z ' X ^ f(z). We stipulate, in addition to X\ being of type H that, furthermore, 
u(z, Afc) = for all Xk{z) for which £(z, A) = 0. Under this additional assumption, we see that in 
fact 

i{X Xx u{z, Xi) + a{z)u{z, A*) = Xf.u(z, A*) (33) 

and we have proven that 
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Theorem 5.2 If X\ is a vector field of type H, u(z, Aj) = and f € C^(D + ), then 

gives an exact reconstruction formula for the density f based on the data I a ,ef of attenuated ray 
transforms of f over the integral curves of Xg . 

6 Conclusions 

The method of complexification presented in the preceding allows for a compact unification of the 
inversion formulae given for ray transforms on both Euclidean space [23J and the Poincare hyperbolic 
disc [3]. Extending the class of vector fields amenable to the aforementioned scheme beyond those of 
type H remains an open problem. Since the analyticity properties of the coefficients of the vector 
fields, ensured by the condition H, were used to justify the holomorphy of the Green's function it 
is unclear how one could alter the method in the absence of such a condition, although the recent 
|21j may yield some insight. With that in mind, there remains the question of finding sufficient (or 
even necessary) conditions on the initial vector field being holomorphic after the complexification used 
above. Real-analyticity is perhaps the simplest necessary condition, but presumably there are much 
more stringent ones. There also remains the question of when u(z, Aj) = 0. 

Lastly, we remark that the only symmetry of the equations occurs when A 6 T, which is not where 
the analysis takes place. In fact, we must break the symmetry in order to arrive at our solution and 
find our minima Aj. Informally, this procedure is analogous to the so-called Higgs mechanism for 
gauge-invariant spontaneous symmetry-breaking of a complex scalar field used in the Standard Model 
of particle physics [121 [2S] ■ 
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